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ABSTRACT 



Using the low energy effective action of the N = 2 supersymmetric SU(2) Yang-Mills 
theory we calculate the free energy, Q, at finite temperature, both in the semiclassical 
region and in the dual monopole/dyon theory. In all regions Q depends on both the 
rS I temperature T and the appropriate moduli parameter, and is thus minimized only for 

specific values of the moduli parameter, in contrast to the T = case where the energy 
vanishes all over the moduli space. Within the validity of perturbation theory, we find 
that the finite temperature Yang-Mills theory is stable only at definite points in the 
moduli space, i.e. for a specific value of the monopole/dyon mass or when the scalar field 
expectation value goes to infinity. 



1. Introduction 

In the last few years, our knowledge about the behavior of supersymmetric non-abelian 
gauge theories has increased dramatically, both for A^ = 1 and N = 2 supersymmetry 
(SUSY) |l], 0]. One of the main achievements has been the derivation of the low energy 
dynamics, expressed in terms of a different, but equivalent, theory. Under certain circum- 
stances this makes it possible to describe the original strongly coupled theory by a dual, 
weakly coupled theory and provides a natural generalization of Dirac's "electric-magnetic 
duality" . 

The classical N = 2 SU(2) gauge theory has a continuos manifold of inequivalent 
ground states. This classical moduli space is parametrized by the gauge invariant quan- 
tity u = Tr(/)^ oc a^, and for a 7^ the original gauge symmetry is broken down to U(l). 
Although there are no quantum corrections that lift the degeneracy among the ground 
states, the quantum theory is still quite different from the classical one. At large ex- 
pectation values the theory is weakly coupled due to asymptotic freedom, and the low 
energy degrees of freedom are just the massless excitations of the microscopic theory. 
Naively, one would expect an enhancement of massless fields at a = 0, where the full 
symmetry is restored. However, in the quantum theory a 7^ all over the moduli space, 
and thus the symmetry group is never enlarged. Instead there exist points where mas- 
sive monopoles/dyons become massless [^. Thus, at different points in moduli space it 
becomes necessary to use different theories to describe the low energy properties. 

Duality in this context is not a symmetry of the theory but rather a way to describe the 
dynamics at a strong electrical coupling in terms of e.g. a weakly interacting, magnetically 
charged theory, where the monopole appears as a fundamental particle. 

Although the N = 2 SU(2) gauge theory is not directly relevant for phenomenological 
studies, it provides a valuable theoretical framework where nontrivial issues can be tested; 
in that respect the supersymmetric theory is worthwhile to study and use as a toy model. 
For example, being an asymptotically free, four-dimensional nonabelian gauge theory, 
several of its fundamental properties, such as confinement via condensation of magnetic 



monopoles p| and chiral symmetry breaking [^, are also found in QCD. 

Another interesting aspect of field theories (including supersymmetric ones), is their 
decisive role in determining the evolution of the early universe, that is their behavior at 
finite temperature. Previous studies of the finite temperature effects in SUSY theories 
have focused on e.g. certain fermion condensates in the non-linear SUSY cr-model P] 
and symmetry non-restoration effects in softly broken SUSY theories at high tempera- 
ture [0. Apart from being interesting in its own right, it is not inconceivable that the 
study of supersymmetric theories at finite temperature could provide additional qualita- 
tive information of a more general nature. Therefore, it is of interest to study also the 
finite temperature effects on supersymmetric theories. More specifically, since the N = 2 
SUSY Yang-Mills theory has a very rich and complicated structure at T = 0, that never- 
theless is rather well understood, one could hope that this theory provides some nontrivial 
information at finite T as well. 

At finite temperature, several properties of a theory may change dramatically from the 
T = case. In particular, a preferred nontrivial minimum of the effective potential can 
be altered as the temperature is raised. In ordinary, non-supersymmetric gauge theories 
this is a well known feature, and in e.g. the electroweak standard model the SU(2) x 
U(l) symmetry is believed to be restored at a critical temperature T^ ^ 200 GeV [0. In 
the N = 2 SU(2) Yang-Mills theory, the existence of a manifold of ground states depends 
crucially on a delicate balance between the perturbative effects from fermions and bosons, 
as expressed in the non-renormalization theorem p. However, at finite temperature 
the bosons obey Bose-Einstein statistics whereas fermions obey Fermi-Dirac statistics. 
Thus, the cancellations due to supersymmetry do not in general apply to the temperature 
dependent parts and one can not a priori expect the finite T effective potential to remain 
unchanged. 

As the SUSY Yang-Mills theory is asymptotically free, one would expect the high 
temperature phase to consist of a weakly interacting gas of "electrical" particles. Another, 
more interesting prospect is the changes in the flatness condition that may occur already 



at low temperatures. A crucial observation regarding the low temperature effects is that 
the ffatness of the potential is sensitive to any disturbance. In contrast to most other 
theories, e.g. the electroweak standard model, this means that the stable points of the 
theory may be altered completely at any temperature T > 0. Thus even a very small 
heating of the system is capable of driving the theory to some specific minima in the 
moduli space. An interesting consequence is that if the temperature initially is very high 
and the "electric" plasma phase is cooled, then eventually the system would be forced 
to one of these minima. Note that since different minima correspond to a completely 
different behavior at low energy , once a minima is chosen the corresponding dynamics is 
fixed and persists even at very low temperatures. This paper is devoted to a study of the 
existence of such low temperature minima. 

To study the low temperature effects on the theory, we will use the low energy effec- 
tive action and calculate the free energy perturbatively. To be able to trust perturbation 
theory the omitted corrections of course have to be small. Imposing this condition we 
will calculate the free energy both in the semiclassical region (large \a\) where the original 
"electric" coupling is small, and in the vicinity of the points at which the monopole/dyon 
becomes massless, where the dual "magnetic" coupling is small. At a given low tempera- 
ture, the free energy is then minimized with respect to the appropriate moduli parameter 
and the stable points can be determined. 

The paper is organized as follows. In the next section we recall some of the main 
results obtained for the low energy N = 2 SUSY SU(2) Yang-Mills theory. Section 3 
is devoted to a calculation of the free energy in the semiclassical region and in section 
4 we analyze the thermal effects in terms of the dual variables, around the monopole 
singularity. Our conclusions and comments on the results are given in section 5, while 
our conventions and some technical details are referred to the appendix. 

2. N = 2 SUSY SU(2) Yang-Mills theory at low energy and T = 

In this section we set up the model and recapitulate the main results at zero temper- 



ature 0. For a comprehensive review, see e.g. p|, |T0|, |TI . 



The N = 2 SUSY SU(2) Yang-Mills theory without matter fields consists of a chiral 
N = 2 field fulfilling a reality condition [l^. In the N = 1 superspace formulation, the 
action can be written as, 

1 



'S'ym = — ImTr / d'^xr 



d^WWa + 2 fd'^d^<l>e'^^<^ 



(2.1) 



where W"' is a vector multiplet containing a vector boson A^ and a Weyl fermion A, and 
the chiral multiplet $ contains a complex scalar (p and a Weyl fermion x- The trace is 
taken over the internal gauge indices, with all fields in the adjoint representation and 
the generators normalized according to Tr(T"T^) = 6°'''/2. The function r is given by a 
combination of the coupling constant g and the theta-angle 6, 

6 4:Tci , , 

^ = 7r + —- 2.2 

271 g'^ 

This theory is asymptotically free, and at the level of perturbation theory the /5-function 
is determined completely by the one-loop contribution P, |l3l , 

(^{g) = -^^i ■ (2.3) 

As in QCD, dimensional transmutation takes place and introduces a scale parameter A 
that governs where the coupling becomes strong. 

Writing the classical action ( p.l| ) in components and eliminating the auxiliary fields, 
the action S'ym contains a potential for the scalar fields, 

r(0t,0) = -LTr([0t,0p) . (2.4) 

This quartic potential is obviously non-negative and is minimized when 0^ and (p com- 
mute. Up to a gauge transformation we can take 

0=2^' (2-5) 

with a an arbitrary complex constant and o"^ the usual Pauli matrix. The gauge invariant 
quantity that labels the set of inequivalent vacua, the moduli space, in the quantum 



theory is given by 

u = (Tr(02)\ . (2.6) 



For a 7^ a super-Higgs phenomenon takes place and the SU(2) symmetry is broken 
down to U(l), whereas the N = 2 supersymmetry is still maintained. 

Now consider the low energy part of the spontaneously broken theory, valid for energies 
E <^ A. When \a\ ^ A, the theory is weakly coupled and the massive particles are very 
heavy compared to the scale E. As a result, to lowest (zeroth) order in a~^ the massive 
particles can be neglected and the low energy effective action simply consists of a free 
massless abelian theory, sitting in a single N = 2 chiral multiplet. Effective interactions 
between the massless particles are suppressed by powers of a and appear when the effects 
of the heavy modes are taken into account. Actually, all terms up to two derivatives or 
at most four fermions can be determined completely (for a proper way of counting in a 
momentum expansion of a supersymmetric theory, see Hl^). In the A^ = 1 formalism, the 
most general U(l) gauge invariant effective action up to two derivatives, and compatible 
with N = 2 supersymmetry is, 
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^cff = Im — / d'^x 
Stt 



d^F"W'^Wa + 2 fd^d'^^F' 



(2.7) 



where F{^) is the so called prepotential and F = dF/d^, F = d'^F/d^'^. The unknown 
function -F($) is then fixed uniquely at the perturbative level by combining the non- 
renormalization theorem |^ and the requirement that the effective Lagrangian reflects 
the 11(1)/? anomaly. Together with the nonperturbative instanton corrections, this gives 

"$2 



' ^.2 



F($) = — <|.^loj 

ZTT 



'A2\2" 



+ ^\E^n(^j ^ (2-8) 



A2 
where Cn are coefficients that can be determined from the exact treatment of the problem 



T3| or by direct instanton calculations K 



Having derived the low energy action in the large \a\ region, one now has to find a 
description involving the correct massless fields that is valid when the coupling becomes 
strong. In the classical theory, the point a = corresponds to an enhancement of the 



symmetry, U(l) -^ SU(2), and all the fields become massless. In the quantum theory, 
a single point where extra massless fields appear is not admissible, because in that case 
a^ would be a global parametrization of the moduli space, which it cannot be 0. Con- 
sequently, there has to be at least two singularities at m 7^ where additional particles 
become massless. Moreover, the appearance of massive vector bosons becoming massless 
at certain points leads to a conformally invariant theory in the infrared limit, which con- 
tradicts u ^ 0. Therefore, Seiberg and Witten 0] conjectured that there should be exactly 
two singularities at ±A^ in the strong coupling region, and that these two points should 
correspond to a magnetic monopole (at +A^) or a dyon (at — A^) becoming massless. 
The dyon has electric and magnetic charge (ne,?T-m) = ^(l) ~1)) while the charge of the 
monopole is {ne,nm) = ±(0, 1). The correct degrees of freedom in the low energy action 
hence include the massless fields that couple locally to the monopole or dyon, together 
with the effective interactions that results from integrating out the light soliton. This 
conjecture leads to a consistent picture and has been justified further later on (see e.g. 



The solitons that become massless at ±A^ exist already at the semiclassical level, since 
the original spontaneously broken SU(2) theory admits soliton solutions. These solitons 
fall into a short multiplet and their mass is given by the BPS-condition, 

M = V2\Z\ , (2.9) 

where Z is the central charge of the supersymmetry algebra. Now, since the relation 
between the central charge and the mass is protected from quantum effects, the mass 
formula ( p.9| ) remains valid even in the full quantum theory. The central charge can be 
calculated explicitly |jTB|, with a resulting mass formula for the BPS states, 

M = V2\ane + aDnm\ ■ (2.10) 

Here a^ = dF/da is the variable dual to a. In other words, the dual theory in the vicinity 
of M = +A^ is a magnetic version of A^ = 2 QED, where the light monopole has a mass 
Mm = \/2\aD\ and couples to the dual photon. When m ^ A^ we have a/) ^ and the 



monopole becomes inassless. The other singularity, where the dyon becomes massless, 
is related to the first by a discrete Z^ symmetry that originates from a nonanomalous 
subgroup of U(l)ij. At the point u = — A^ the dyons are massless, corresponding to 
(a — a^) -^ 0. 

Integrating out the monopoles, the resulting effective theory is reminiscent of ( p.7|) , 
but in terms of the dual fields Wj^ and $/). Given the /3-function for N = 2 SUSY QED, 

1 



Pd = ^9l , (2.11^ 



the dual prepotential becomes 



FB($B) = --^$^log 

An 



(1)2 

A2 



+ A^E^nf^) , (2-12) 



n=l ^A 



where the coefficients c^ can be calculated from the exact solution \^\. 



3. The free energy in the semiclassical region 

In this section we will calculate the free energy for large \a\, based on the low energy 
effective action. 

As has already been noted, for |a| ^ A the effective action at energy scales E <^ A 
contains only the remaining massless fields. These fields are not charged with respect 
to the unbroken U(l) and are thus free to first order. All possible interactions between 
the massless fields can be expanded in local, effective vertices and the interactions are 
suppressed by powers of the masses of the heavy particles, which are proportional to 
\a\. Also, the coupling constant is small due to asymptotic freedom which suppresses the 
higher order terms further, albeit only logarithmically. 

If the temperature is low, i.e. T ^ \a\, the probability of having a real massive 
particle in the heat bath is suppressed by the Boltzmann factor e"''*'/"^ <^ 1. Hence 
we can safely neglect the thermally excited massive particles, and need only to take into 
account the exchange of virtual, massive particles between the massless and otherwise free 
ones. Therefore, the dominating contributions to the free energy come from the massless 

8 



particles and their effective interactions. Furthermore, when the temperature is low we 
do not expect any large logarithms involving the temperature, so g = g{\a\) <^ 1. 



The component form of the effective action can be obtained from equation (|2.7| ) in a 
straightforward manner. Eliminating the auxiliary fields and including all vertices that 
contribute to the free energy up to |a|^^, the action is given by. 



Si^ =ld'x\ (d.p^Wp) - '-x^x -\m- If^'F,. + 



167r2 



p[ _P 
a^ a 



F'"'F,,+ 



/lU 



+ 



47r2 



^3 



p[^P 
a^ a 



j'tJo/A' 



{d.p^nd'-p) - '-xpx - \m - If'-'f,^ 



;i + 75)X(M + ^^7'^(1 + 75)A(a^p) - V2 AS'^^(1 + 75)%^^-' 



+ 



9 



IGTT^at 

,4 






+ 



^X7'^(l + 75)X(5/.P^) + *A7'^(1 + 75)A(9^pi-) _ ^2 xS^'^(l - 75)AF^. 
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647r2 

,6 



9 

'1671 



- {X(l + 75)X A(l + 75)A} + -^ {x(1 - 75)X A(l - 75)A} 
(^— ) {A(1 + 75)A x(l + 75)X - X(l + 75)A A(l + 75)x} + 
+ (^-^) {A(1 - 75)A x(l - 75)X - X(l - 75)A A(l - 75)x} + 
+ (^){A(1+75)A A(1-75)A + X(1 + 75)XX(1-75)X + 
+ 2x(l + 75)AA(l-75)x} 



(3.1) 

where the fermions A, x have been written as four-component Majorana spinors, p is a 
complex scalar and F/j^^, is the field strength for the vector boson, with F^^, the dual of 
Ff^i^. To arrive at the action ( |3.1| ) all fields have been rescaled by the common factor g. 
Finally, the instanton contributions have been neglected, since they are further suppressed 
by powers of (A/|a|) <ti 1. 

Since the free energy is a static quantity we will use the imaginary time formulation 
and apply by now well established techniques for such finite temperature calculations 
P^ , ppj] . As is well known, at finite temperature T, the Euchdean space is restricted to 
xG3fJ'^,0<X4</9 = T^^ and the bosonic (fermionic) fields are periodic (antiperiodic) 
under X4 — *> X4 + /?. 

To first order, we have a gas of free, massless particles and the corresponding contri- 



bution to the free energy density Q is just the blackbody radiation, 



fi(i) 



TT 



2rp4: 



12 



(3.2) 



At two loops, there are several terms of order |a|~^ in equation ( |3.1| ) that could con- 
tribute. However, it is easily seen that apart from the vertices proportional to either 
{p^/a^) or {p/a), the order |a|~^ loops from all the other terms vanish trivially, since the 
remaining integrands are either odd or contain an odd number of 7-matrices. Concerning 
the contribution from the diagrams involving {p^a^) and (p/a), they are all of the form 
depicted in Fig. 1. Since the momentum integrals are ultra-violet divergent, they have 
to be regularized. Employing a regularization method that respects the supersymmetry, 
such as dimensional reduction^ [|2ll , |2^ , and using the method outlined in ||23| to factorize 
the integrals, the result of the derivative interactions can be written in the form. 



fi(2) 



oc 



(27r)' 



± 



(27r)'^p2 



+ 



d\ 
(2^ 
d^ 



loo+e 



dko 



loo+e 
ioo+e 



2TTi e^^^o =F 1 
dpo f l\ 



X 



(27r)^ j-ioo+e 27ri \p'^ J e- 



/3po 



(3.3) 



We have inserted an arbitrary mass parameter p to get the correct dimensions, and the 
sum over the Matsubara frequencies Un-, ojm have been written as contour integrals, with 
fco = iujn, Po = ioJm and (3 = T~^. The second bracket in the free energy (|3.3|) is due 
to the scalar field, where uj^ = 2imTT and in the first bracket the upper sign holds for 
bosons and the lower for fermions, where Un = (2ra + l)7rT. Now, in the first bracket the 
first term vanishes due to the definition of integration in dimensional regularization. In 
the second term inside the first bracket we can close the contour in the right half plane 
of the complex /cQ-plane, and since the integrand has no poles there is no contribution 
from this term either. As a consequence, the expression in equation ( p.3| ) vanishes and 
n^'^^ = 0. Thus there are no r^/|ap corrections to the free energy. 



^The conclusions below are also true when the dimensional regularization method is used. 
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scalar field propagator 

scalar, fermion or 
vector field propagator 



Fig. 1. Nontrivial diagrams in n^'^\ at order 



a\ 



-2 



Expanding the action (|2.7] ) further to include higher orders in |a| ^, it is easy to see 
that there are no loop corrections to the blackbody radiation that are suppressed by a^. 



Therefore, the leading correction term to the free part in ( |3.2| ) has to be at least of order 
0(T^/|a|^). These terms will also be of order g^ or higher in the coupling constant. 

However, the order T^/|a|'^-terms do not only originate from an expansion of the effec- 
tive action in ( |2.7| ). Apart from this holomorphic piece, there are also non-holomorphic 
contributions that start at the level of four-derivative terms P3, ^^, BH]. In the A^ = 1 



superspace formulation the non-holomorphic part of the effective action is | |14|| , 

+ if$$ll{41^°W^aT^aW-8iy"D„$W^D>t}] , (3.4) 

where we have defined 

A- = log(i)log(*), (3.5) 

and the subscripts on K denote derivatives with respect to that argument. Rewriting 
equation (|3.4| ) in component form and taking into account only the terms up to |a|~^, we 
find that all possible divergent terms vanish one by one when the dimensional reduction, 
or dimensional regularization, procedure is used. The remaining terms in the effective 
Lagrangian that contribute to the finite temperature effects are given by, in the Weyl 
two-component notation. 



L 



(2) _ -9^ 



"^ 167r2|a|4 



/ p^lv rp X 2 

lid,p^)id^p){d,p^){d''p) + 2p^p{d,d,p^){d>^d''p) + 2 [-^] + 
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+ 2 ( ^^ 1 + {2{x^^d,x){d.X^''x) + 2(x5mX)(x^^x) + (x ^ A)} + 

+ {z(a'^x^'^x)(V)(^-p) -^(x^^5'^x)(5.P^)(M +^P(^/^^^P^)(X^'^^^X)- 
- ip\d,d,p){d^x^''x) + ^(x^^5.x)i^,5F'^' - *(9,x^^x)i^M5i^'^' + (X - A)} - 

+ (Aa'^a^9,x)(x5^A) - {xdA){\a^a>^d,x) ] + O ((^Vkl') , (3-6) 

where all the fields have been rescaled with a factor g, in accordance with conventions 
used for the holomorphic action (|3.1|) . The terms of order g^ appear when the lowest 
order equations of motion for the auxiliary fields, derived from the holomorphic action 
( p.7|) , are substituted back into the action (|3.6| ). Since the lowest order terms in the 
non-holomorphic part are of order g'^, it is clear that the leading correction terms to the 



free-field result come from equation ( |3.6| ), and not from the holomorphic part (|3.1| ). 

The contributions to the free energy come from two-loop diagrams of the form shown 
in Fig. 2, where each loop, independently of the other, is any of the bosons (p, A^) or 
fermions (A, x)- 




Fig. 2. Two-loop contribution in D.^'^\ at order \a\ ^. 



Evaluating the different diagrams, we find the following contributions. 



fermion-fermion terms: n|?> = 1?^ 

^ 43200a^ 



scalar-scalar terms: fiif-* = — -- , 

2700|a|4 ' 

vector- vector terms: \1}J = —- — — -r , 

2700 a^ ' 



scalar-fermion terms: ^2 
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(3)^77rVT 
""^ 5400|a| 



• vector-fermion terms: ^y = , 

• vector-scalar terms: fl^^J = — -r . (3.7) 

1350|ap 

As expected, since the physical degrees of freedom for the complex scalar field and the 
vector boson are the same, their contributions are identical. Also, the fermions and bosons 
obey different statistics, and the result of having a fermion loop or a boson loop should 
not be equal. However, as is apparent from the terms in (|3.7| ), they all give results that 
add up instead of tending to cancel each other, as is the case in SUSY theories at T = 0. 
This additivity of the thermal effects is in accordance with earlier observations on the 
finite T behavior of SUSY theories [^ . 

Combining the free-field part in ( |3.2| ) and adding all the terms from ( ^.7] ), we finally 
have 

12 192|a|4 \^ |a|4 J ^ ' 

Since we have found a dependence on both T and \a\ in Q, we can now regard the 
temperature as being a fixed external parameter and minimize the free energy with respect 
to \a\. As is clear from equation ( p. 81) , the free energy decreases as |a| becomes larger, and 
since a larger \a\ corresponds to a smaller coupling constant, the approximations made 
in (|3.8|) become even better. Thus, the prediction from the perturbative calculation is 
reliable and the free energy tends towards its minimum value when \a\ -^ oo. This run- 
away solution in the semiclassical region is somewhat reminiscent of the superpotential 
solution in A^ = 1 SU(NJ Yang-Mills theories with N/ (Nj < NJ flavors |28|, although 



the origin of the behavior is of course completely different. It is a rather interesting result 
that the semiclassical region is well defined at T = 0, but tends to an unphysical solution 
for any T > 0. 

However, it should be stressed that the above conclusion is only trustworthy to the 
extent that higher order corrections in g can be neglected, and is strictly speaking only 
true in the far asymptotic region. In fact, if the order g^ contribution from (|3.6|) and the 
order g^ correction from the holomorphic part (|3.1|) are non-zero, and if the numerical 
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coefficient in front is large, one can not exclude that the non-asymptotic region has 
another behavior, with an additional stable point at large, but finite \a\. 

4. The free energy in the dual theory 

In this section we will use the dual formulation of the theory and calculate the free 
energy around the monopole singularity. 

As |a| decreases, the theory moves away from the semiclassical region and the system 
becomes strongly coupled. Although the condition \a\ ^ T might be fulfilled, the strong 
coupling renders a perturbative calculation in the original "electric" variables meaningless. 
Fortunately, the powerful tool of duality provides a mean to perform the calculation in 
certain regions of the strong coupling regime. When entering the strong coupling region 
the lightest BPS-states consist of the ±(0, 1) monopoles around m = A^ and the ±(1, —1) 
dyons for u ~ — A^. 

As stated earlier, around the point u = A"^ the low energy theory is a magnetic version 
of A^ = 2 QED, where the magnetic monopoles have a mass Mm = v2\aEi\ and become 
massless as -u — * A^. The dual coupling constant qd is small as long as the characteristic 
energy scale E is small, E <^ A. 

Restricting the calculation to T <^ \a\ and T <^ A still permits an arbitrary relation 
between the temperature and the monopole mass. Two extreme limits can be calculated 
analytically, namely T ^ Mm and T ^ Mm- In the first limit monopoles easily become 
thermally excited, whereas in the latter one has to use an effective action where the 
monopoles have been integrated out. 

Now, consider first the case when {aol <C T (and T ^ \a\, T ^ A). Under these 
conditions we are in the weak coupling region, and the heat bath consists of both the 
massless N = 2 chiral multiplet (where the "dual photon" is the vector boson) and the 
massive but light N = 2 monopole hypermultiplet. At this temperature the contributions 
from any other particles, such the dyons, are neglectable since their mass is proportional 
to \a\. Thus we are in a high temperature phase with respect to the light degrees of 
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freedom. 

Using A^ = 1 superspace, the Lagrangian of the theory is, 

1 






57r 



-Im <^ To 



<fd W^Woa + 2 / d^rff <I>D$ 



'D 



+ Qoe^''^QD 



£e 



V2Qd^dQd + mQoQ 



D 



+ h.c. \ , {4.1] 



where m = \/2aD- Wd, $d belong to the N = 2 chiral multiplet and Qd, Qd to the 
hypermultiplet, where Qd and Qd have opposite gauge quantum numbers. After some 
rescahngs, the component form of the Lagrangian ([4.1|) becomes, 



.(D) 
^eff 



1 



--F.^F^'' - -X^X - -xH - ^i^m^ + ( V)(^'^P) + pA.0)n^^0) + 



\^iTM^A 



4 



+ [b^^)\m) + ^ U(l - 75)A0 - A(l + 75)V^0t _ A(l - 75)^0+ 



+ V(l + 75)A0t 



9d_ 

V2 



pip{i + 75)?/^ + pV(i - 75)^^ + 0x(i + 75)^+ 



+ 0V(1 - 75)X + 0^(1 + 75)X + 0M1 - 75)^' 



+ 



+ V2gD fM^(0^0 + 0V)(P^ + P) - V2(7dpV(0V + 0^0) 



9h 

2 



-J. ~\2 



(0V + 0t0) +M^V^z^-M^(0V + 0t0) 



(4.2) 



where -^ = {d^ —igoA^^)!^, F)^ = d^ —igoA^ and D^ = d^ + igDA^^. The Lagrangian has 
been written in a four-component form, and (p, A^, x. A) belong to the chiral multiplet 
whereas {(j),(f),ip) come from the hypermultiplet (and ip is a Dirac spinor). 

As before, the lowest order contribution to the free energy is just the blackbody ra- 
diation. Denoting the free energy density in the dual theory by Qd at T :^ Mm, the 
contribution from the free part is. 



Q 



(1) 

D 



2t^4 



n'T 



+ 



m 



MlT 



+ 0(M;^log(M^/7rT) 



(4.3) 



6 4 37r 

where the free energy for the massive particles have been expanded in powers of {Mm/T). 
As is clear from equation ( |4.3| ), the noninteracting part of the free energy becomes smaller 
when Mm -^ 0. Based only on the part of the free energy displayed in ( |4.3| ), one would 
conclude that Mm = represents a minimum. However, as will be clear, the interaction 
effects are not insignificant. 
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Taking into account also the interactions from the Lagrangian 



the next contri- 



butions to Qn are two-loop effects, which are suppressed by 0{gjj). Using the method 



developed in ||29| , |30| , their contributions can be calculated in a power series in (Mm/T). 
We find, 



fermion-vector terms: 
(2) _59lT' i, , 



n 



18 /M^ 



, fv 



288 [ 
scalar-vector terms: 



57r2 \ T2 



7 



1 (M„ 



fi 



(2) 
Dsv 



144 \ 
Yukawa terms: 

o(2) _ 25^j,T^ I 

^y " 288 1 

scalar-scalar terms: 

79lT' i 



1- 



24M„ 



12 /M^ 



7Tn 



77r2 \ r2 



VttT 



7 - - + log 



[attTJ 



^ 48M^ ^ 



72 /M2 



25T7r 257r2 \ r2 



37 
4 



3 9 , 

— I log 

8 12 ^ 



ttT 



+ 77 log(4) 



n 



(2) 
Dss 



144 



1- 



30Mm 

7Ttt 



M2 



77r2 \ T2 



57 



23 

y 



Slog 



UttT/ 



(4.4) 



where 7 is Euler's constant and terms of order OlgjjM^T) and higher have been omitted. 
Adding all the different parts in equation ( [4. 41 ) together with the free-field result in 
we have, 



2t^4 



n 



D- 



'~6 



29^ 

247r2 



2 rp2 



13gj,M^T'_ ^ MiT 



, 1 , /M„ 

H log --— 

16 ^ V ttT 



+ ^log(4) 



247r ■ 4 



l + M(l + ^+ 
7r2 132 16 



37r 



+ o(^|,m^t,m; 



(4.5) 



The total free energy in ( [4. 51 ) displays some remarkable features. When T ^ M^ 
the relevant scale is set by the temperature, so we expect that g^j^ = g'jj{T). Now, for a 
sufficiently small Mm/T ratio, the term proportional to g]^MmT^ is comparable to M^T'^. 
In other words, the interactions may change the minimum (at Mm = 0) found from the 
blackbody radiation only. From Q^ in ( [4.5|) , the dominant mass dependent contribution 
QMm in ^D is 



Q 



Mr, 



4 247r / "" 



(4.6) 
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and therefore the position of the extremum of the free energy is, 

M.„^MI. ,4.T) 

In fact, this point represents the minimum. Hence the free energy is minimized when 
equation ( [4.7| ) is satisfied, that is when the zero temperature mass is of the same magnitude 
as the thermally induced "electric" mass. By the Z2 symmetry this implies a similar result 
to be valid in the region where the dyon is the lightest massive particle. 

We now turn to the opposite extreme in the dual region, i.e. when M^ ^ T . We 
also require Mm <^ A to be able to trust the perturbation theory. When the parameters 
fulfill the above conditions the thermally excited monopoles only contribute with a factor 
g-A/m/T Q£ Qourse, the exponentially small term reflects the probability of having real 
monopoles in the heat bath. As in the semiclassical region, it is clear that the dominating 
effect from the monopoles is due to the virtual exchange of monopoles between the massless 
particles, and it is appropriate to use a low energy effective action. 

The holomorphic part of the effective action is given by equation (|2.7|) , but in terms 
of the dual variables IVd, $£> and with the dual prepotential ( |2.12| ). The contribution 
from the blackbody radiation to the free energy f2^ is, 

^ff = --^ . (4.8) 

As in the previous section, we have to include the interactions to incorporate the mass 
parameter, in this case IodI- In the semiclassical region the holomorphic part did not 
contribute to the leading correction terms in the free energy, and this result persists in 
the dual theory as well. The reason is simply that the expansion of the dual prepotential is 
completely analogous to the case in the semiclassical region, with the expansion parameter 
a in the semiclassical region replaced by a" -^ —'ia'^, except in the coupling constant. 
This in turn is a consequence of the fact that the coefficient in the /3-function of A^ = 2 
QED differs by a factor —1/2 from the SU(2) case, as a comparison between equation 
(^ and (Pl| ) shows. 

However, just as in the case of ordinary, non-supersymmetric QED one should expect 
higher order interactions coming from an Euler-Heisenberg type of effective Lagrangian 

17 



3l|| . In fact, such a supersymmetric extension was discussed in [^, and corresponds 
to nothing but the non-holomorphic contribution to the effective action. In that sense, 
the non-holomorphic action in the semiclassical region is the supersymmetric extension 
of a nonabelian Euler-Heisenberg action, whereas for the magnetic monopoles it is an 
extension of the abelian case. 

As in the holomorphic part, the effective action in the dual region is given by replacing 



a^ -^ —2a]y in ( p.6[ ) (see also |2^). Therefore, the lowest order correction to the free 
energy, when Mm ^ T, is 

The total free energy in the dual region, valid for Mm ^ T is then obtained by combining 
equation ( ]i7^ ) and (^^), 

^D = ^^^ — ^ + f^ . 4.10 

12 384|aB|^ Vl^sl / 

Minimizing the free energy we thus find that the theory prefers a smaller mass, a statement 

that of course is valid only as long as both Mm ^ T and the omitted higher order 

corrections are neglectable. Again this result is transferred to the region around the dyon 

singularity by the Z2 symmetry. 

To summarize this section, we have found that when Mm ^ T the free energy becomes 
smaller as the mass decreases. When the mass is sufficiently small, the theory enters a 
high temperature region where T ^ Mm- As expected, the mass continues to decrease, 
but rather unexpectedly only down to Mm oc g^T . If the mass initially is smaller, it will 
increase instead, up to the minimum of the free energy at Mm oc g\)T . Note though, that 
as in the semiclassical region we can not exclude another minimum in a region where higher 
orders in the coupling constant become relevant. Since the dual coupling increases as the 
mass becomes larger, such a minimum could be present when the theory is sufficiently far 
away from -u = A^, but still in the regime where the dual theory is well defined. 

Using the Z2-symmetry, we also conclude that the same situation as described above 
takes place around the singularity at m = — A^, where the dyons are light. 



5. Conclusions 

We have calculated the free energy in the low energy N = 2 supersymmetric SU(2) 
Yang-Mills theory, using the low energy effective action description, and would now like 
to make some comments and remarks about our results. It should be emphasized that 
some of the conclusions are naturally rather speculative at the present stage. 

In the semiclassical region where \a\ ^ A, we find that the lowest order correction 
terms come from the non-holomorphic part of the effective action, and as a result the free 
energy decreases as the moduli parameter \a\ increases. The stable point corresponds to 
the limit \a\ —>■ oo. Using the dual description in the strong coupling regime, we find that 
the theory flows towards a stable point where the monopole mass Mm is of the order of 
the thermally induced mass, M^ oc (^^T. The discrete Z2-symmetry relates the monopole 
singularity to the region where the dyons become the light solitons, and therefore implies 
that the monopole results are valid for the dyons as well. Thus the scenario depicted in 
Fig. 3 emerges, where region III corresponds to the semiclassical region \a\ ^ A, and the 
two other regions correspond to a neighborhood where the dual theory is valid; region II 
to where the soliton mass M^ satisfies Mg ^ T and region I to T ^ Mg. 



in II -y-.i 



/f '- \ 



^:X' " ^ 



X '■' X 



Fig. 3. Directions where the free energy decreases, at different parts in moduh space. 

Dotted hnes denote the stable minima. 

As mentioned earlier, one can not rule out the existence of other minima, positioned 
at large but finite \a\ and at a large value of the soliton mass, Ms, in the dual theory. In 
order to establish whether such minima actually appear, one would have to calculate the 
logarithmic corrections coming from higher orders in the coupling constant, but having 
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the same suppression of powers of mass as the dominant effect (i.e. \a\ ^ or M^ ^). At 
this stage the presence of such minima can neither be proved, nor excluded. 

In the asymptotic region, where the minimum is reached for \a\ -^ oo, the theory does 
not seem to be well defined in the finite temperature case. Nevertheless, the large \a\ 
region could still be meaningful in e.g. a cosmological scenario. This would be attributed 
to the fact that ultimately higher energy scales must become important and hence may 
provide a stabilization of the theory, as has been speculated to be the case in the A^ = 1 
theories with a similar runaway solution at T = ||28|| . 

In the strong coupling region the stable points correspond to |a/)| oc g'j;)T, so that the 
phase oi ao is irrelevant. The position of the minima at {aol is not likely to be modified by 
higher order terms in the low energy effective theory. However, the possible implications 
of the heavy particles, whose mass is proportional to \a\, remain unclear. For instance, if 
the monopoles are strictly massless the value of the free energy is higher than at the point 
Mm = (13(y'|,T/127r). But since the slope at the origin is zero there is no tendency to float 
away from the massless case, and the theory seems to remain at a local maximum. The 
ultimate fate of the massless theory must then be determined by the effects of the heavy 
particles. 

Finally, we would like to comment on the different phases that exist at low tempera- 
tures. If the temperature initially is much higher than any of the internal scales A and u 
in the microscopic theory, we expect to have a nearly free gas of the "electrical" particles. 
As the temperature is lowered, a description of the theory in terms of the low energy 
effective action becomes valid. Whatever the nature is of this transition, the theory will 
end up in one of the minima that exist at low temperatures. The mechanism that drives 
the theory to a specific minimum presumably depends on how the scalar field expectation 
value develops as the temperature is lowered. This is certainly a complicated process at 
an energy scale where standard methods such as perturbation theory is not applicable. 
Without reliable information it is obviously very dangerous to make any predictions, but 
it is not unreasonable to expect that different regions in space could fall into different min- 
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ima. In that case, domain walls would form and remain until the temperature vanishes. 
Needless to say, this scenario is presently by no means rigorous. 
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Appendix 

In this appendix we fix our notation and illustrate the calculations of the free energy 
with a simple example. 

Starting in Minkowski space-time, our metric has the signature 

<7'^^ = diag(+ ) , (A.l) 

and the antisymmetric tensor has eoi23 = +1- The 7-matrices satisfy {7^^, 7''} = 2g'^'' 
and are given by 

where a^ = (1 , a), a^ = (1 , —a) and a are the usual Pauli matrices. Also, 

S^^ = I [r, 7l • (A.3) 



From the actions ( pTl|) and (|4.2| ) we find the following propagators. 



fermions: 



scalars: 



1^ — 771 ' 

i 
9 9 ' 



vector bosons: — , (A. 4) 

pi 



where m is the possible particle mass. 
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To calculate the free energy, we first write the logarithm of the partition function Z 
in Minkowski space-time and then make the following replacements, 

2nS{po) -> -tl36^„ , (A.5) 

where (3 = T~^, ujn = 2miT for bosons and uo^ = {2n + l)7rT for fermions. Labelling the 
expression that follows when the Minkowski rules and the replacements in equation ( [A.5| ) 
have been used by {—ipVJ^), with V denoting the volume of the system, the free energy 
density is, 

logZ 



n 



(3V 



iT . 



(A.6) 



To give a simple illustration of the procedure outlined above, we will calculate the 
contribution from the interaction term 



-(D) 



9l 



a. 



in the dual theory. We find. 



^E 



Sp 



a. 



1 



(A.7) 



(27r)3 (^u,^Y-p-Ml 
We now write the sum over the Matsubara frequencies as a contour integral |]19|, pO|, 



(A.8) 



Tj2f{po = iuJn) 



-ioo 2Tli 



f{po)+fi-Po) 



± 



ioo+e 



ioo+e 



dpo 

27ri 



fiPo) + f{-Po) 



oPPO 



Tl 



(A.9) 



where the upper sign holds for bosons. Omitting the temperature independent piece, we 

have 

r 1 2 

2/"'^^ 1/ 1 \ 

= '^[JJ2^^ ^p + Ml leVp-+M^ - l)\ ' ^ °^ 

by closing the contour in the right half-plane and picking up the pole at po = Jp"^ + M^. 
When T ^ M^ we can expand the remaining integral in [Mm/T) by use of the following 
result. 



dxx"^ 



TC 



y^ , 2/^ fi 



\^x'^ + y'^ \f,^/x^+y^ 



1. 



T-f-^ti2-^-'°Hi^)f-^o(^''-(*-") 



y 
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Equation ([A.10|) then finally becomes, 



n ^- 
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